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A procedure based on the semiclassical approximation for high energy levels is developed to 
yield solutions to the classical equation of charge motion and to the Bargmann-Michel-Telegdi spin 
equation. To this end, exact solutions to the Klein-Gordon and the Dirac-Pauli equations are used. 
The essence of the procedure under review is that the quantum state of a charged particle in a 
homogeneous magnetic field is represented as a superposition of states corresponding to the neigh- 
boring energy levels. As a consequence, the behavior of the expectation values of the momentum 
and spin operators with respect to the resulting nonstationary wave function (packet) strictly obey 
the classical equations of charge motion and spin precession. 
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I. INTRODUCTION 



In the semiclassical approach, the charge motion is described in a classical way, and quantum transitions are 
calculated on the base of matrix elements [Q. This method is known to be fairly in the ultrarelativistic case at 
ultrahigh electron energies. For example, at energies of ~ 2.5 GeV, which are typical for current storage rings with a 
characteristic parameter Hp ^ lO'' Oe cm the electron energy levels amount to n ~ 10^''. For radiative transitions to 
lower energy levels, An/n ^ 10^^. In this case, according to the uncertainty principle, all quantum processes occur 
in a region as small as a few angstroms (see Q for details). Clearly, this is a negligible small value for an orbital 
radius of about 10 m. Hence, the classical, or more precisely, semiclassical concept of the electron trajectory can be 
used here. 

In the theory of synchrotron radiation, the semiclassical method allows all quantum corrections to the synchrotron 
radiation power, including the effects with spin-flip to be calculated This method is straightforward to use. 

What is more, it provides a visual picture of physical processes occurring at ultrahigh electron energy (recoil effects, 
radiative spin self-polarization, spin magnetic-moment radiation, mixed charge and magneton radiation, radiation 
associated with the anomalous part of electron magnetic moment, etc.). 

Despite of the impressive success of the semiclassical synchrotron-radiation theory rigorous mathematical substan- 
tiation of this method based on exact solutions to the Klein-Gordon equation or to the Dirac-Pauli equation is, to 
our knowledge, lacking heretofore. The basic statements of the theory were, in fact, postulated on the basis of the 
uncertainty principle (see 

In this paper, solutions to the classical equation of charge motion and to the Bargmann-Michel-Telegdi (BMT) spin 
precession equation Q are constructed in the framework of quantum mechanics, using nonstationary wave functions 
(see general statements in |^). The later are, in their turn, exact solutions to the Klein-Gordon equation or to the 
Dirac-Pauli equation for a charged particle moving in a homogeneous magnetic field. A simplified version of this 
method was used in |Q to reveal a physical pattern underlying the behavior of longitudinal electron-spin polarization 
in a magnetic field. This approach is adopted in ||^ to describe quantum neutron spin-flip transitions. 

Our task is to explore the time evolution of the expectation values of the momentum and spin operators for a 
charged particle moving in a homogeneous magnetic field. To demonstrate the capability of the method discussed, let 
us examine the problem at hand first. In Sec. ||, we will consider a scalar particle and compare the time evolution of 
the expectation value of the momentum operator for this particle with the solution to the classical equation of charge 



motion. These are found to differ by a factor. Sec. HI will focus on a spin-i particle. In addition to evolution of 



the expectation value of the momentum operator, we will discuss the behavior of the expectation value of the spin 
operator. The latter also differs from the solution to the BMT equation by factor only. Finally, in Sec. IV, we will 
introduce a procedure that will help to eliminate the difference between the quantum and classical approaches. To 
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this end, the state of a particle will be represented by a wave packet so that the time evolution of expectation values 
of the corresponding operators with respect to this nonstationary wave function will coincide with that predicted 
classical mechanics. 



II. MOTION OF A SCALAR PARTICLE 

Let us assume that state of scalar particle is formed as a wave packet by superposing exact solutions to the Klein- 
Gordon equation corresponding to three neighboring energy levels n — 1, n, and n + 1: 

fi 



*(r,t)= ^ A^V,„(r)exp(-^TOoc2B™i), 

m—n—l 



where V'm(r) are the stationary solutions to the Klein-Gordon equation in a homogeneous magnetic field (see Appendix 
0) and Am are the expansion coefficients. 

To determine the coefficients A.„i in the semiclassical approximation (n ^ f), the probability that a particle will be 
found at each of the levels is assumed to be the same. For n 3> 1 this is natural assumption. Normalizing the wave 
hmction 5'(r,t) to unit probability yields a relation of the form 



n+l 

E 



A^Am — 1. 



In what follows an explicit form of Am is immaterial but different relations between these coefficients are of great 
importance. In particular, we need the following equation 

SIkG = A*^n-l + K+lAn = I 

Now let us calculate the expectation value of the momentum operator P in the state ^(r, i). In the first step 
we calculate matrix elements of (i/'m' IPIV"™) type with respect to stationary solutions to the Klein-Gordon equation. 
Then, using the semiclassical condition {n ^ 1) for each matrix element [see Eqs. ([A^)] we write (^|P|^')j in explicit 
form to give 

i 2 
{Px)t — T:'mocb± [SIkg exp(ia;i) — SI^g •3xp(— iwt)] = —-mocb±smujt, 

1 2 
~ -TOoc6j_ [21kg exp{iujt) — SI^q exp{—iujt)] — -mocb± cosujt, 



Here 



Z — \t>n^\ — tin 



is the frequency which, in essence, coincides with the cyclotron frequency of rotation of a classical particle in a plane 
perpendicular to the magnetic field vector. In this approximation, _B„ —^ 7, where 7 is the Lorentz factor. 

This leads us to conclude that the expectation values differ from the corresponding classical solutions to the equation 
of motion of a charged particle only by the factor (2/3) in the {Px)-t and {Py)^ components. Note that the same result 
was also obtained in [^ for an ordinary one-dimensional quantum harmonic oscillator. In the discussion below we will 
show that this factor can be made as large as unity. 
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III. CHARGE MOTION AND SPIN PRECESSION OF A DIRAC-PAULI PARTICLE 



When considering the Dirac-Pauh particle the spin polarization of the particle must be taken into account. A 
superposition of exact solutions to the Dirac-Pauli equation corresponding to three neighboring energy levels n — 1, 
n, and n + 1 has the form 

n+l 
^ m—n—l 

where 'i/'m(r) are the stationary solutions to the Dirac-Pauli equation (see Appendix ^ and Ai^m are expansion 
coefficients. Note that in this case, allowance is made for polarization states defined by a quantum number ^ = il- 
Let us assume that a given longitudinal-polarization state at time t = 00 for an electron located at level m is of the 
form 

(cr • P) ^ A^m^mC = KnAQmlpmC- (2) 
C 

Having solved Eq. (0) we find the eigenvalue 



In the semiclassical approximation, it can be supposed that for a ll va lues of m, the spin coefficients are the same, i.e.. 



Ci{mX) — CiinX) and, moreover, Am ~ A„ [see Eqs. ( |B3[) and (|B4|)]. Then we have 



Let us also require that the wave function ^'(r, t) be normalized to unit probability. As a result, we get 

ji+i 

J2 E Al^Arn = l- (3) 
m—n—l 

In what follows we need the relations 

" 2 
E Km^i™+i+^-i™^-i™+i)= 3, (4a) 

m—n—l 

n 

m—n — l 

" 2 k 

= E ^llm^lm+l = o— T> (4b) 
m—n—l 

n+l 

2l3= E ^U^-lm-^^, (4C) 

hi -hi 

7n— n — 1 

n+l ^2 

2t4= E (^Im^lm - ^-lm^~l™) - ^VVT- (^^^ 
m—n—l 



^It is known that, with regard to the anomalous magnetic moment of an electron, the operator (cr ■ P) is not an integral 
of motion. 
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To simplify the calculations we suppose that, as an electron makes a transition from one level to another, the 
projection of the momentum of the electron onto the direction of the magnetic field is conserved. In addition, random 
deviations of the orbital center of the electron due to radiation are neglected. These restrictions imply that b'^ = bz 
and s' = s (see Appendix 

Now, by analogy to a spinless particle, we find the expectation value of the operator P in Dirac-Pauli particle 
state (see the corresponding matrix elements in Appendix As a result of straightforward but cumbersome 
calculations, we derive, in view of these restrictions, the following equations: 

i 2 
{Px)t — 7:''TT'0cb± exp{iujt) — 21* cxp(— iwt)] — — -moc5_L sinwt, (5a) 

1 2 

{Py)f = -mocb± exp(iwf) — 21* exp(— icji)] = -mocbj^ costot, (5b) 

{Pz)t = niocbz. (5c) 

Here 



CO 



is the frequency which in the BMT approximation (the charge motion is independent of the spin precession) coincides 
with the cyclotron frequency of rotation of an electron in the plane perpendicular to the magnetic-field vector. In this 
approximation, i?„^ 7. 

A similar procedure applies to calculations of the expectation value of the spin operator S"^ (see 0) to yield 



Taking into account Eqs. (Q), we obtain 



{Sx)t = — (±{cosujtsmQat + b sin Lut cos D, at), (6a) 



2 

{Sy)^ = C± (sin ti't sin r^ai — b cos ujt cos il, at), (6b) 

3 

iSz)t = ^C. + ^-4^C±cosnat, (6c) 



(^°), = jCz + S„cya COSilat. (6d) 

Here Cz = 1 ~ C±i C-L = 2k{k^ + 1)~^ are constants which, as we will be seen later, are analogous to some constants 
in the classical spin theory. Thus, as in the classical approach, the total precession is determined not only by the 
cyclotron frequency w but also by the anomalous frequency 

TOqC^ fJ^aH b 

"a — 7 — [BnC — Bn^i) — 2 



2'K moc 



Note that in the BMT approximation, the term comprising the anomalous magnetic moment in the spin operator 
is immaterial for result (o). 



We see that by analogy with the spinless charged particle, the results obtained in Eqs. (pa\), (5b), (|6a|) and (6b) are 
slightly different from the corresponding solutions to the classical equations of charge motion and spin precession by 
the factor (2/3). 
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IV. SEMICLASSICAL CORRESPONDENCE PRINCIPLE 



To construct a more precise semiclassical theory of c harg e motion and spin precession, we will consider a wave 
packet involving exact solutions to Dirac-Pauli equation (B2) corresponding to closely spaced energy levels^. In this 



case, wave function (]1|) can be written in the following form: 



^'(r, i) = X! X! ^C"iV'mc(r)exp(--^TOoc^B„^t), (7) 

where 3? is a set of values of the principle quantum number. Let us assume that each value of m in Eq. (Q) belonging 
to the set is much higher than unity, i.e., to ^ 1. At the same time, we have 

n > iV = TOmax - TOmin > 1, 

where TOmax and TOmin are maximum and minimum values belonging to 3?, n is the quantum number from 3? defining 
some energy level to which a given classical trajectory corresponds. Let us further assume that the other energy levels 
from the set are symmetric about level n. 

The next natural step in our method is to derive the coefficients The central idea of the derivation is the 

same as in Eq. (||), i.e., 

^Im = KA_ijn, y ' y ^ j4^„jj4^„j — 1. 

C meSR 

Then, relations (Q) can be easily extended to the case under consideration, namely, 

rrimax-l N -1 

2li= = — r^, (8a) 



lm+1 



m— mj-nin 



E.^. N — 1 K, , , , 

= (8b) 

2I3 - E ^ImA-im = (8c) 

mGSR 

2t4 = E (^^m^l™ - ^-l™^-l™) = 4-T- (8d) 

Following the above procedure of calculating the expectation values of the momentum and spin operators and using 
Eqs. (^) we obtain 

A^- 1 

\Px)t = j^rnQcb^ smujt, 



N - 1 



(Py)j ~ — — — mocb^ cos Lot, 



^It can be shown that there is a simpler form of this construction for the motion of a spinless charged particle. 
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N - 1 

{Sx)t = — (^±{cos cut sin flat + b sin ujt cos flat), 

N -1 

{Sy)^ — — — C-L (sin sin fiat — b cos cut cos ^ at), 

\S^)t = 5 + — ^C±cosrJai, 

Thus, for iV ^ 1, the difference of {Px)ti {Py)tJ {^x)t, and {Sy)^ from the results obtained by classical theory is 
eliminated, and the time evolution of the expectation values (P)^ and {S'^)^ are made to coincide with the solutions 
to the corresponding classical equations |p^ . 

It is easy to see that the anomalous magnetic moment of an electron has effect not only on the time behavior of 
spin projection onto the direction of motion (longitudinal polarization) |^ but also on the behavior of spin in the 
plane perpendicular to the magnetic-field vector, i.e., on the time evolution of {Sx)^ and {Sy)^. 

V. CONCLUDING REMARKS 



It should be noted that derived in Sec. IV satisfy the relations given by the classical spin theory: 

Moreover, the procedure under review can be extended to the tensor operator |^ 

* = —P3{o- X P) + -^pi(<T X H), 

moc moc^ 

fl = <T + —P2{(T X P) + -%H, 

which, in the BMT approximation, is related to by the formula 

TOqC 

Thus, the expectation values {S^^)^ and (P)j as well as (II^'')^ calculated following the procedure discussed in 



Sec. IV obey the classical equation of charge motion and spin precession. The approach introduced in this work 
substantiates the use of the semiclassical method for gaining an insight into physical phenomena associated with 
high-energy charged particles. In particular, the semiclassical approach can be the basis for a classical model for 
certain purely quantum processes. For example, it can be used to establish a relationship between spin-flip transitions 
and classical spin precession . 
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APPENDIX A: 



Motion of a scalar particle in a constant and homogeneous magnetic field H — (0, 0, H) is given by the Klein-Gordon 
equation 

where E = ihd/dt is the energy operator, P = p — (e/c)A is the kinetic momentum operator, A = {—Hy/2, Hx/2, 0), 
e = —Co < is the charge, and toq is the rest mass of the particle. 

It is known that in this case, the wave function ip{r, t) is to be the eigenfunction of the energy operator E as well 
as z-components of the kinetic momentum operator P and orbital angular momentum operator L = r x p onto the 
direction of magnetic field vector, i.e., 

Elp = En^, Pz^ = Pzip, Lzip = hltp. 

Then, in terms of cylindrical coordinates (r, ip, z), which are the most suitable ones in this case, the solution to the 
Klein-Gordon equation has the form [0 

^(r,t) = ^^^ e-'^"*/''e'P--/'^e''^/„..(p). (Al) 



The function In,s{p) in Eq. (Al) is defined by the Laguerre polynomials Qiip) with help of the following relation 



Here n = I -\- s = 0,1,2,... is the principle quantum number, s = 0, 1, 2, . . . and / = 0, ±1, ±2, . . . are radial and 
azimuthal quantum numbers. 

The energy of a particle is written as 



E,, = moc'Br, = moc'Jl + bl + bl, bi_ = 2,/ -^(n + i), 
V V mgc^ 2 

where — Pz/moc is the projection of the momentmn onto the direction of the field and fiQ = eQh/2mQC is the Bohr 
magneton. 

Matrix elements of the momentum operator of scalar particle have the form 

{i^m'\Px\lp7n) = ■^mocb±{6rn'-m, + l " ^m'^m,-l), (A2a) 

(V'm'lA/IV'm) = ■^■mocb±{S„i'^,n,+l + <5m'-m,-l), (A2b) 
i'ipni'lPzl'ipTn) = mocbzSrn'rn- (A2c) 



APPENDIX B: 

In the paper we use the solutions to the Dirac-Pauli equation 

= H = c{a-P)+ pamoc^ -I- PaPsi^- ■ H), 

ot 

where P = p — (e/c)A is the kinetic momentum operator, A = {--Hy/2, Hx/2, 0), H is the external magnetic field, 
pa = (Q;/27r)/io is the anomalous part of electron magnetic moment, e = — eo < and nio are the charge and rest 
mass of electron. 
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It is known that in this case complete set of commuting operators characterizing the quantum state of the 
particle consists of the Hamiltonian H, the projections of momentum operator P and total angular momentum 
operator J = L + ■^hcr onto the direction of magnetic field vector which, for definiteness, is oriented along the Z-axis, 
and polarization operator. According to |Q, in a homogeneous magnetic field the operator 



P2{(T X P 3 + 

mgC c 



is taken as spin integral of motion characterizing the polarization of electron spin relative to the magnetic field vector. 
Then we have 



(Bla) 



n,^ = C = ±1- 



(Bib) 



In terms of cylindrical coordinates (r, ip, z), the wave function being a solution to the Dirac-Pauli equation and 
satisfying Eqs. (Bl) has the form |Q 



where 



^(r t) — _^L£2^_e"* 



( Ciln-l,s{p) ^ 

I iCdnApy"^ ) 



(B2) 



Here n = I + s = 0, 1, 2, . . . is the principle quantum number, s — 0,1,2,... is the radial quantmn number, and 
I = 0, ±1, ±2, ... is the azimuthal quantum number. In dimensionless form, the spin coefficients Ci are determined as 



BnC 



'1 



nC 



BnC 



where the energy of electron is given by the expression 

EnC 



Bn. 



< 



1 + 6i + C 



PaH 

niQC^ 



(B3) 



(B4) 



bz = , Ox = 2W -n, 6 = 

moc V "^oc 

Let us consider exact forms of all necessary matrix elements of the momentum and spin operators. 

z 

{'4'm'C'\Px\'<Pmc) = ■^'mocb±{Sm'-m, + l " Sm' -m,,-l)5c' ( , 
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(V'm'f'l-PjIV'mC) ='moCbz5m'mScC^ 



ii'm'cl^xl^mc) = ■^{i' ~ C(<5m'-m, + l " <5m'-m,-l))(<5m'-m,+l — (^m'-m, -l)<5-C'C' 



{ipm'cl^vli'mc} = ~ C(^m'-m,+l - <^m'-m,-l))(^m'-m,+l + ^m'-m,-l)<5-C'C> 



[1] J. D. Jackson, Rev. Mod. Phys. 48, 417 (1976). 

[2] Synchrotron Radiation Theory And Its Development, edited by V. A. Bordovitsyn (World Scientific, Singapore, 1999). 
[3] V. A. Bordovitsyn, V. S. Gushchina, and A. N. Myagkii, Nucl. Instr. Methods A 405, 256 (1998). 

[4] V. A. Bordovitsyn, I. M. Ternov, and V. G. Bagrov, Usp. Fiz. Nauk 165 1083 (1995) [Sov. Phys.-Usp. 38, 1037 (1995)]. 

[5] V. Bargmann, L. Michel, and V. L. Telegdi, Phys. Rev. Lett. 2, 435 (1959). 
[6] D. M. Fradkin and R. H. Good, Rev. Mod. Phys. 33, 343 (1961). 

[7] I. M. Ternov, V. G. Bagrov, R. A. Rzaev, and Yu. I. Khrrionko, Izv. Vyssh. Uchebn. Zaved. Fiz. 6, 111 (1964). 
[8] G. G. Cabrera and M. Kiwi, Phys. Rev. A 36, 2995 (1987). 

[9] I. M. Ternov, V. G. Bagrov, and V. A. Bordovitsyn, Izv. Vyssh. Uchebn. Zaved. Fiz. 4, 41 (1967). 
[10] V. A. Bordovitsyn, Izv. Vyssh. Uchebn. Zaved. Fiz. 36 11, 39 (1993) [Sov. Phys. J. 36, 1046 (1993)]. 
[11] A. A. Sokolov and I. M. Ternov, Radiation From Relativistic Electrons (American Inst, of Physics, New York, 1986). 



9 



